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We present a study of the effects of dispersion on the properties of photon pairs generated by 
type-II, collinear spontaneous parametric downconversion. Specifically, we take into consideration 
the effects of a chirped pump, as well as of dispersive propagation of the photon pairs. We present 
expressions for the joint amplitude both in the spectral and temporal domains, as well as for the 
chronocyclic Wigner function of heralded single photons, which fully characterizes the single photon 
spectral (temporal) properties. On the one hand, we show that unwanted effects of pump chirp in 
terms of the heralded single photon duration can be suppressed for states designed to be factorable 
and spectrally elongated. On the other hand, we show that pump chirp constitutes an effective tool 
for the control of the degree of photon-pair entanglement. We show that when frequency-entangled 
photon pairs propagate through a dispersive medium, entanglement can "migrate" between the 
modulus and phase of the joint temporal amplitude. 

PACS numbers: 42.50.Dv, 03.67.Bg 



I. INTRODUCTION 

The ability to generate two-photon states with spe- 
cific properties is essential for quantum information pro- 
cessing applications. Specifically, it has been determined 
that factorability is required for heralding of pure single 
photons [1], a crucial ingredient for quantum computation 
with linear optics [2 . In previous work, various recipes 
have been proposed, some of which have been demon- 
strated experimentally, for the generation of factorable 
photon pairs through spontaneous parametric processes 
IlHlEElIZllHlElllOlIIT]. In particular, in a recent 
experiment, we have demonstrated the effectiveness of 
asymmetric group velocity matching for the preparation 
of factorable two-photon states, and hence of heralded 
single photons with a high degree of purity [U E]. In 
this experiment we have exploited the concept of engi- 
neering the photon pairs at the source so as to avoid the 
need for filtering, thus achieving a high flux of factorable 
photon pairs. This scheme and similar ones have one as- 
pect in common: they require a broadband pump, often 
one in the form of a train of femtosecond-duration pulses. 
Indeed, a monochromatic pump precludes factorability, 
leading to a state with a high degree of spectral entan- 
glement. As is also the case with classical optics, the 
use of femtosecond-duration pulses, greatly increases the 
vulnerability to dispersive effects. Naturally, dispersion 
leads to a significant effect on the resulting photon-pair 
properties. It is crucial to understand this effect for the 
design and implementation of photon pair sources tai- 
lored for quantum information processing applications. 
In this paper we present a comprehensive study of the 
effects of dispersion, in particular of second order disper- 
sion, on the properties of photon pairs generated by type- 
II, collinear parametric downconversion. We include in 
our treatment the effects of both: a quadratically chirped 



pump and propagation of the signal and idler photons 
through a medium with second-order dispersion. 

Previous work has analyzed some aspects of dispersive 
effects on photon pairs generated by parametric down- 
conversion. Thus, for spectrally anti-correlated pho- 
ton pairs it has been demonstrated that: i)the correla- 
tion time remains unaffected if the signal and idler pho- 
tons propagate through separate dispersive media char- 
acterized by opposite-signed quadratic dispersion [13], 
and ii) the shape of the Hong-Ou-Mandel interference 
dip is immune to quadratic dispersion experienced by 
one of the generated photons [M]. The latter effect has 
been proposed as the basis for a quantum optical co- 
herence tomography scheme which benefits from disper- 
sion cancellation [15]. The temporal broadening of pho- 
ton pairs due to dispersion which underpins these effects 
has been been directly observed [16 . Recently, non-local 
temporal shaping of the signal photon through chirping 
of the idler photon has been demonstrated [TTl [TS] . The 
above results rely on a monochromatic pump, and there- 
fore do not consider or exploit chirp in the pump beam. 
In Ref. [19 , it was shown that two-photon states with 
positive spectral correlations, produced by a broadband 
pump, can be immune to dispersion, to all orders. 

In the present paper we analyze the interplay of an 
ultrashort, quadratically chirped pump with quadratic 
dispersion experienced by the signal and idler photons 
generated by the process of type-II, collinear PDC. Our 
motivation for this work is to analyze the various dis- 
persive effects which may be relevant in the design of 
specific sources of photon pairs for quantum informa- 
tion processing applications. We derive expressions for 
the resulting joint amplitude, both in the spectral and 
temporal domains, explicitly including the effects of dis- 
persion. On the one hand, we show that the negative 
influence of pump chirp in terms of the attainable her- 
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laded single photon temporal duration can be reduced, 
or even eliminated, through a dispersion suppression ef- 
fect which occurs for a state designed to be factorable 
and spectrally-elongated. On the other hand, we show 
that pump chirp can in fact be a powerful tool for con- 
trolling entanglement in photon pairs. We study the re- 
lationships between: i)the degree of entanglement, ii)the 
heralded single-photon purity and iii)the expected Hong- 
Ou-Mandel interference between two photons generated 
in different crystals. We study the spectral (temporal) 
properties of heralded single photons. Finally, we dis- 
cuss an effect which occurs during propagation of the 
PDC photons through a dispersive medium, by which 
entanglement can "migrate" between the modulus and 
the phase of the joint amplitude. 



describes a spectral filter acting on each of the signal and 
idler modes. 

The phase matching function can be shown to be given 

by 

^{uji^ujs) = smc[LAk{uJi,uJs)/2]ex.p[iLAk{uJi,uJs)/'2], 

(3) 

in terms of the phasemismatch Ak{uJi^uJs) = kp{uJs + 
^i) — ks{ujs) — ki{(jOi) and the crystal length L. In this 
work we rely both on numerical calculations taking into 
account dispersion to all orders, as well as on analyti- 
cal expressions based on a Taylor series description of 
the phasemismatch. In the latter case, we can write the 
phasemismatch as a function of the frequency detunings 
= — (with ji = s^i) as 



II. THE TWO PHOTON STATE FOR TYPE-II 
PDC UNDER THE EFFECTS OF DISPERSION 

In this paper we study the spectral (temporal) proper- 
ties of photon pairs produced by spontaneous parametric 
downconversion (PDC) pumped by a train of ultrashort 
pulses in cases where the pump, signal and idler fields ex- 
perience dispersion. In particular, we study the effects on 
the two-photon state of: i) quadratic chirp in the pump 
pulses (e.g. due to a dispersive element prior to the crys- 
tal), and ii) propagation of the two photon states after ex- 
iting the crystal through a medium with quadratic chirp, 
such as a fiber. Throughout this paper, we refer to dis- 
persion introduced by optical elements prior and follow- 
ing the crystal, i.e. excluding dispersion introduced by 
the crystal itself, as external dispersion. We will concen- 
trate our discussion on collinear, type-II (with the signal 
and idler photons orthogonally polarized) and frequency 
degenerate parametric downconversion. 

Following a standard perturbative approach, the two- 
photon state produced can be written as 



+ biuf + bpUs^i (4) 

where perfect phasematching is attained at cj^ = cc;^ = cc;c, 
Ak^^^ is the constant term in the series (which we assume 
to vanish), and in terms of the following definitions 

= LW^{2uj,)-kl{u;,)] (5) 

b, = ^[k;{2u,)-k';,{u;,)] (6) 
bp = Lk';{2uj,). (7) 

Here, (with fi = s^i) represents a group velocity mis- 
match coefficient between the pump and the signal/idler 
fields, while bp^s,i represent group velocity dispersion co- 
efficients involving the pump, signal and idler frequencies. 

The broadband pump is described by the spectral enve- 
lope function a{Lj)^ which we model as a Gaussian func- 
tion with bandwidth a 



|ilf) = |vac) +7/ J dujs j dujif{uji,ujs)\u;s)s\^i)i, (1) 

where we have assumed that the two-photon state is 
spatially-filtered so that only those k vectors which are 
parallel to the pump direction of propagation are re- 
tained. Here, |cc;)^ = a]^(cj)|vac) with fi = i^s where 
|vac) is the vacuum, /(coi^cos) represents the joint spec- 
tral amplitude (JSA) and 77 is a constant related to the 
conversion effiency. In the presence of dispersion in ad- 
dition to spectral filters on the paths of the signal and 
idler modes, the JSA can be written as 

(2) 

where ^{uOi^uOs) denotes the phasematching function 
(PMF), a{uJi-\-uJs) represents the pump spectral envelope 
function (PEF), D{u;i,u;s) describes a phase term associ- 
ated with external dispersion, and F/j^^cj) (with ji = i^s) 



Oi{^s + ^i) = exp [-{us + z^O^/cr^] • (8) 

Likewise, function F^iuo) (with ja = s,i) is modelled 
as a Gaussian function with width ap, i.e. F/^^u) = 
exp(— z/^/(7|.). By neglecting cubic and higher order dis- 
persive terms, and also neglecting linear dispersive terms 
(which result in a temporal shift without otherwise modi- 
fying the two-photon properties) , we can express function 
D{vi,Vs) as 

D{vi, Vs) = exp [i(3p{iyi + i/^)^] exp [ipiiyf] exp [ipsJ^'^] 

(9) 

in terms of the group velocity dispersion (GVD) param- 
eters (with /i = p^s^i). If these (3 coefficients are due 
to propagation through a dispersive medium of length 
and characterized by wave number a^^((x;) (placed prior 
to the crystal in the case of the pump and following 
the crystal in the case of the generated photons), then 
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P^i = where " denotes a second frequency deriva- 

tive, evaluated at ujc for PDC light and at 2uJc for the 
pump. Throughout this paper we employ a temporal 
description of the two-photon state, based on the joint 
temporal amplitude (JTA) function /(t^,^^), obtained as 
the two-dimensional Fourier transform of the JSA (see 
Eq. [2|. In terms of the JTA we may write down the 
two-photon state as: 



1^) = |vac)+7? j dts j dUf{U,ts)\ts)s\U)i (10) 

where \t) ^ = a]^(t)|vac), defined in terms of the time- 
domain annihilation operators 



(11) 



the joint amplitude, both in the spectral and temporal 
domains, based on a Gaussian approximation for the sine 
function in the phasematching function. Clearly, this ap- 
proximation ignores the sidelobes associated with the sine 
function and leads to slightly different spectral properties 
even near perfect phasematching, compared to the full 
calculation. However, the resulting joint amplitude ex- 
pressions capture the essential physics of the two-photon 
state and can be exploited for the derivation of analytic 
expressions for the following: i) Hong-Ou-Mandel inter- 
ferogram, relying on two independently-generated her- 
alded single photons (Secc. [TTT| ), ii) Chronocyclic Wigner 
function for heralded single photons (Secc. [V|) , and iii) 
Schmidt decomposition of th e modulus of the joint tem- 
poral amplitude (Secc. VII). Thus, by approximating 
sinc(a::) ~ exp(— 70;^) with 7 ^ 0.193 (this value is se- 
lected so that the two functions have identical full widths 
at half maximum), we can write down the JSA as 



In order to calculate the JTA analytically we use 
the power series description of the phasemismatch (see 
Eq.|4|, up to first order terms. We also neglect the phase 
term in the phasematching function (see Eq. [3|, which is 
justified in many of the cases to be analyzed in this paper 
for which external dispersion dominates over dispersion 
introduced by the crystal. The resulting expression is as 
follows 



f{ti,ts)= iVtexp 
1 



4Co 



X l^erf 
— erf 



1 



(Ci(ti,t,) + 2Co) 



(Ci(ti,t,)-2Co) 



(12) 



Here, 

T^^ (with X,ii = s, i) are similar to parameters Tf,,, with 



,) = iV^exp ^{Ttivt + T^y, + 2TtiV,Vs)\ ■ (15) 
Nuj is a normalization constant, and parameters 

~2 



an additional term in their real parts 



-'Am,/ 



rTi2 



(16) 



The corresponding expression for the JTA is as follows 



fiti, ts) = iVtexp [Qlt^ + ^Itl - 2^litiU] , (17) 

where Nt is a normalization constant and where we have 
defined the parameters with ji = i^s and Vtsi as 



where Nt is a normalization constant, and which is ex- 
pressed in terms of the following definitions 



Co 



I6(f4-T;2f2j 
t^jfjTs - fin) + ujfln - f^Ts 



r< (f f ) — ss'^i '^-^ si'^t'^a ^ -List's 



(13) 



Here, we have defined parameters T^^ (with X^ji = i) 
in terms of their real and imaginary parts 



T:2 



(14) 



where 6x^ represents a Kronecker delta. A similar ex- 
pression, for the case without external dispersion, has 
been presented in [2^. We also present expressions for 



4(T4 -T2.T2J 



(18) 
(19) 



If the joint amplitude is normalized so that 
/ dvs j dvi\ f{vi,Vs)? = 1, |/(^i,^s)P represents a joint 
probability distribution for the emission of photon pairs 
with frequency detunings and Vi. We refer to 
|/(^i,^s)P as the joint spectral intensity, or the joint 
spectrum. It is convenient to define, for later use, a cor- 
relation coefficient S = —agJ{as(Ji) where represents 
the covariance associated with this probability distribu- 
tion and where cr^ (with ji = s^i) represents the standard 
deviation for each of the two marginal distributions. The 
correlation coefficient is constrained to take values within 
the range — 1 < S < 1; in terms of the parameters defined 
in Eq. [16) the correlation coefficient is given by 



^ si,R 



(20) 
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An important class of two-photon states is composed 
of those which are unentangled, or factorable. For these 
states, functions S{u;) and /(cj) exist such that the JSA 
may be expressed as /{coi^cos) = S{u;i)I{u;s)- Factorable 
states are required for heralding of pure single photons. 
In previous work, it has been shown that in cases where a 
symmetric group velocity matching condition is fulfilled 
[3, 21 , symmetric and factorable two-photon states, with 
a round joint spectrum, are possible. Likewise, it has 
been shown that in cases where an asymmetric group 
velocity matching condition is fulfilled, factorable states 
with an elongated joint spectrum are possible. 

For the asymmetric group velocity matching case 
where = k'^ or = (i.e. where the pump and 
signal photons have identical group velocities), the gen- 
eral expression for the JTI (see Eq. 12 ) may be shown to 
reduce to 



\hu.ts)\' 



= Na exp 1^- 
xrect (-^ 



2(l+/32a4)^ 

O 1 ^.H t^7. 



(21) 



where Na is a normalization constant. See the appendix 
for a definition of the rect function and for more details 
on how the general form of the joint temporal amplitude 
(see Eq. 12) simplifies for particular cases. 



Throughout this paper we will refer, as an example 
of interest, to factorable and spectrally elongated two- 
photon states obtained by asymmetric group velocity 
matching (GVM). Another, related, category of state is 
that of factorable, symmetric two-photon states obtained 
by symmetric GVM. In Fig. [l] we present, for illustration 
and future reference, the joint spectral intensity as well 
as the joint temporal intensity for specific examples of 
these two categories of source, in the absence of exter- 
nal dispersion; later figures explore the effect of external 
dispersion. These plots are derived from numerical cal- 
culations, taking into account the full two-photon state 
with dispersion to all orders (see Eq.|2|. In the case of the 
factorable, asymmetric state we have assumed a pump in 
the form of an ultrashort pulse train centered at 415nm 
with a FWHM bandwidth of 5nm (which corresponds to 
CF = 4.65 X lO^^rad s~^). Likewise, we have assumed a 
2cm-long potassium dihydrogen phosphate (KDP) crys- 
tal with a cut angle of 0pm = 67.8°, chosen for type-II, 
collinear, degenerate phasematching. In the case of the 
factorable, symmetric state we have assumed a pump in 
the form of an ultrashort pulse train centered at 757nm 
with a FWHM bandwidth of 15nm (which corresponds 
to cr = 4.19 X lO^^rad s~^). Likewise, we have assumed 
a 2.29mm-long P barium borate (BBO) crystal with a 
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FIG. 1: (a) Joint spectral intensity for a nearly- factorable 
two-photon pair source, which fulfils asymmetric group veloc- 
ity matching (based on a KDP crystal), (b) Corresponding 
joint temporal intensity, (c) Joint spectral intensity for a 
nearly-factorable two-photon pair source, which fulfils sym- 
metric group velocity matching (based on a BBO crystal), 
(d) Corresponding joint temporal intensity. 



III. ENTANGLEMENT AND SINGLE PHOTON 
PURITY 



One way to characterize the spectral (temporal) prop- 
erties of photon pairs is through their interference prop- 
erties. Two photons interfere in a Hong-Ou-Mandel in- 
terferometer (HOMI) with perfect visibility if they are 
indistinguishable. In the case where the two photons 
come from the same source [see Fig[2|a)], indistinguisha- 
bility leads to the requirement that the joint amplitude 
be symmetric. In the present analysis, where we con- 
centrate on the spectral degree of freedom, this means 
specifically that /{coi^cos) = /(cj^, cj^). In the case where 
the two interfering single photons are heralded from in- 
dependent sources, indistinguishability leads to the re- 
quirement that the joint amplitude be factorable. Con- 
cretely, we assume the experimental situation depicted in 
Fig. [2|b), where two single heralded photons each in the 
signal mode of two identical sources, are made to inter- 
fere at a beam splitter. We assume that the joint spectral 
amplitude functions for the two sources are given by the 
functions fi{uji^ujs) and /2(^i,^s)- In our analysis, let 
us begin with the expression for the four-fold coincidence 
rate as a function of the delay between the two interfering 
photons in the Hong-Ou-Mandel inteferometer ^22j 
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Rc{r) = l— / dcji / dLJ2 / dcos / <^^4/i (^i, ^2) 
Jo Jo Jo Jo 



(22) 
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FIG. 2: (Color online) Schematic diagrams for the single- 
source [panel (a)] and double-source [panel (b)] Hong-Ou- 
Mandel interferometer setups. While these diagrams show, 
for graphical clarity, PDC sources emitting into non-collinear 
signal and idler modes, all calculations in our paper refer to 
type-II collinear sources where the two modes can be distin- 
guished by polarization rather than by optical path. 

Here, we have assumed that each of the joint amplitude 
functions is normalized. The subscripts 1-4 in the above 
expression refer to each of the similary-labelled paths in 
Fig. [2] Assuming that the minimum coincidence rate (i.e. 
the lowest point in the interference dip) occurs at zero 
delay, the interference visibility V, defined as the depth 
of the dip normalized by the background coincidence rate, 
is given by 



= dcj did' Psi(ll},Ll}')Ps2{^' 

Jo Jo 

= Tr(p,ip,2). (23) 

Here, we have defined the partial density operator pg^ 
for the interfering photon from each of the two sources 
(/i = 1, 2), obtained by tracing over idler frequencies; the 
corresponding matrix elements are 

P5^(cj,cj') = {ij\Tri{p^)\ij') 

duJof^{uJo,uj)f^{uJo,uj') (24) 



/ 

^0 



in terms of the two-photon density operators p^ = 
and the partial trace over idler frequencies Tr^. 
If we assume that the two sources are in fact identical 
then, from Eq. 23 , V = Tr(p^) which represents the pu- 
rity of the heralded, interfering single photons. There- 
fore, from an experimental measurement of the two- 
source Hong-Ou-Mandel interference dip, we may infer 
the purity of the interfering single photons. 

Let us now express the joint amplitudes (assumed to 
be the same for both sources) in terms of a Schmidt de- 
composition [23] 



(25) 



where are the Schmidt eigenvalues and Un{uJi) and 
Vn{uJs) are the Schmidt functions. Note, that the Schmidt 
decomposition can also be performed in the temporal do- 
main, in terms of the corresponding eigenvalues An ^ and 
temporal Schmidt functions Un\ti) and Vn\ts) 



(26) 



It is straightforward to show that the reduced density 
matrix elements can be expressed in terms of the spectral 
Schmidt functions as 



(27) 



which in general represents a mixed state. Calculating 
the visibility from Eq. [23j assuming that both interfering 
single heralded photons are identical {psi = Ps2 = p), we 
can furthermore show that 



y = Tr(p2) = ^A: 



2 _ 



K- 



(28) 



Here we have exploited the orthogonality of the 
Schmidt functions Un{oo). As may be appreciated from 
Eq. [28j the HOMI visibility in the specific case where the 
two interfering photons are identical corresponds to re- 
ciprocal of the Schmidt number K. (The Schmidt num- 
ber K can be used to quantify the degree of entangle- 
ment; while K = 1 represents a factorable, i.e. unentan- 
gled state, K 1 represents a highly entangled state). 
Therefore, according to Eq. |28j from an experimental 
measurement of the two-source Hong-Ou-Mandel inter- 
ference dip, we may likewise infer the degree of entan- 
glement of the two photon states. Let us note that in 
the derivation of Eq. [28] we have assumed an ideal situ- 
ation for which possible experimental imperfections can 
be disregarded. In a realistic situation, Hong-Ou-Mandel 
interference will be influenced by many experimental fac- 
tors, with the implication that it will no longer be possi- 
ble to infer the degree of entanglement from the HOMI 
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visibility. In this case, however, a lower bound for the 
single-photon purity or an upper bound for the degree 
of entanglement (rather than the actual values) can be 
inferred from the Hong-Ou-Mandel measurement. 

From the above discussion it is clear that the gener- 
ation of pure, heralded single photons requires a source 
factorable photon pairs, i.e. one for which there is a sin- 



gle term in the sum of Eq. 25 (and therefore, also in the 
sum of Eq.[27l). 



For the two photon state expressed in terms of the 
phasemismatch up to first order and in terms of the Gaus- 
sian approximation (see Eq. 15), the four- fold coincidence 
rate for the two-source HOMI arrangement is given by 



i?c(r) = l-Vexp(-rVAr2) 
in terms of the visibility V 



V = 



and the dip temporal width Ar 



rp2 rp2 

-^ii,R-^ss,R 



(29) 



(30) 



At' 



rp2 rp2 

ii,R ss. 



rr2 
-^ii,R 



(31) 



Note that according to Eq. [28) there is a reciprocal 
relationship between V and K so that the Schmidt num- 
ber may be determined as = 1/V from Eq. 
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It may 

be seen that the visibility reaches unity, as expected, if 
the mixed term responsible for spectral entanglement T^^ 
vanishes. Note that while the dip shape (visibility and 
width) exhibits a dependence on the pump chirp parame- 
ter, it is independent of the dispersion experienced by the 
signal and idler photons after exiting the crystal. This is 
as expected, since the dip shape is a signature of the type 
and degree of entanglement present, which cannot change 
through lossless propagation of the signal and idler pho- 
tons after exiting the crystal. 

We illustrate this behavior in Fig. [3) in the context 
of a factorable, spectrally elongated source based on a 
KDP crystal (with the same parameters as assumed for 
Fig. [ija)). In particular, we show the expected HOMI 
dip for four different levels of pump chirp (see figure cap- 
tion). While the dips shown in Fig. [3]^a) were calcu- 
lated numerically from the full two-photon state, those 
in Fig. ^h) were calculated through the analytic expres- 
sion in Eqns . [29| through [31] In general terms, the greater 
structure in the joint spectrum, e.g. derived from the sine 
function, in the full calculation results in lower visibil- 
ities compared to the approximate analytic calculation. 
As expected, a larger degree of pump chirp leads to lower 
HOMI visibilities, signalling a greater degree of entangle- 
ment in the two-photon states. In particular, no pump 
chirp leads to perfect visibility (according to the approx- 
imate calculation) and to the highest attainable visibility 
{V = 94%) in the case of the full calculation. 



^0.5 



"^0.5 



(arm 






r 

I: 



-10 





(ps) 



10 



FIG. 3: (Color online) Fourth- fold coincidences (normalized 
to a unit background) vs temporal delay between two sep- 
arate heralded single photons interfering through a Hong- 
Ou-Mandel interferometer. These curves result from: (a) 
a numerical simulation, based on Eq. |22| and (b) the ana- 
lytic expression in Eq. |29| The different curves in each panel 
correspond to the following levels of pump chirp: i)pp = 
(shown in red), ii) f3p — 1.19 x 10~^^s^ (shown in green), iii) 
/3p = 2.39 X 10-^^^ (shown in blue), iv) /3p = 3.58 x 10"^^^ 
(shown in orange) and v) [3p — 4.77 x 10~^^s^ (shown in ma- 
genta); in this family of curves, lower visibilities correspond 
to higher levels of pump chirp. 



In the next section, we further explore the important 
relationship between pump chirp and the resulting degree 
of entanglement. 



IV. CONTROL OF ENTANGLEMENT: USE OF 
A CHIRPED PUMP 

The information pertaining to dispersion in the two- 
photon state, as may be seen from Eqns. |2] and [9j is con- 
tained in phase terms. Note that the phase introduced by 
propagation of the signal and idler photons through dis- 
persive media is factorable into separate signal and idler 
contributions and therefore does not contribute to en- 
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tanglement. In contrast, the phase term associated with 
a chirped pump (e.g. due to transmission of the pump 
beam through a dispersive medium before reaching the 
crystal), in general cannot be factored into signal and 
idler factors and therefore does contribute to entangle- 
ment. This contribution of pump chirp to entanglement, 
in the spectral domain is necessarily in the form of phase 
entanglement [24 , and has no effect on the joint spectral 
intensity. However, phase entanglement in the spectral 
domain can, as we will discuss, translate into expHcit 
correlations in the joint temporal intensity. Of course, 
apart from external dispersion considered in this paper, 
other experimental parameters such as crystal length and 
pump bandwidth can have a profound effect on the re- 
sulting degree of entanglement; these effects have been 
studied elsewhere (see for example Refs. [3 and [20 ). 

The phase term due to pump chirp, linear in the sum of 
frequencies a;^ + cjg , in general manifests itself as an elon- 
gation of the joint temporal intensity, along the direction 
given by U + ts, when plotted in the times of emission 
space {ti^ts}. We illustrate the effect of pump chirp, 
presenting as an example the case of nearly-factorable 
two-photon states (in the absence of pump chirp) ob- 
tained through asymmetric group velocity matching. In 
Fig. [4j we present the joint temporal intensity expected 
for a KDP crystal, with the same parameters as as- 
sumed for Fig. [ija). In this case, the pump group ve- 
locity equals that of the signal (ordinary- wave) photon. 
Fig. [4]^a) shows a plot of the Gaussian term in Eq. 21 
vs the times of emission ti and tg, assuming /3p = 0. 
Fig. ^h) shows a corresponding plot of the rect function 
in Eq. |2l] Fig. shows the joint temporal intensity 
|/(^i, ^s)P in the case of no pump chirp, obtained from the 
product of the functions plotted in the previous two pan- 
els. From Eq.[2l] the effect of pump chirp becomes, clear: 
the width of the Gaussian function in Fig.|4]^b) increases, 
thus revealing more of the diagonal structure provided by 
the rect function in the joint temporal intensity. This be- 
havior is clear from Fig. ^d) , where we have plotted the 
JTI for a pump chirp parameter Pp = 3.04 x 10~^^s^ 
(corresponding to propagation through a 6cm thickness 
of fused silica). 

The phase term in the phasematching function (see 
Eq. |3|, which is not taken into account in the analytic 
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and 



15), in general 



le two-photon 



expressions of section |n| (i.e. Eqns. 
can contribute to the entanglement in i. 
state. Indeed, the mixed term proportional to Us^^i in 
Eq.[4]leads to a phase term in Eq. 15 which cannot be fac- 
tored into signal and idler contributions. This small effect 
can be graphically appreciated by comparing Fig. [ijb), 
which takes into account numerically the crystal phase 
term in the JSA, with Fig. |4jc) calculated analytically 
from Eq. |2l] which does not take into account the crys- 
tal phase term in the JSA. We can see that portions of 
the JTI involving higher values t^, emitted near the sec- 
ond crystal face, involve a broader signal photon, relative 
to portions of the JTI involving lower values of ti. This is 
a consequence of temporal elongation of the pump pulse 




-1 1 
ti (PS) 



FIG. 4: This figure shows the effect of pump chirp on the 
temporal properties of the two-photon state. In particular we 
show, plotted as a function of the signal and idler times of 
emission: (a) the Gaussian function in Eq. 12 for the case 
of no pump chirp, (b) the rect function in Eq. [12] (c) the 
square modulus of the product of the functions in the previous 
two panels, (d) same as in (c), but for a chirped pump with 
/3p = 3.04 X 10"^^s^ 



(which leads to a corresponding signal mode temporal 
elongation) as it propagates through the crystal. This 
intra-crystal dispersion and its effect on the resulting de- 
gree of entanglement, may be compensated for by pump 
chirp, involving dispersion opposite in sign to that ex- 
perienced by the pump beam in the non-linear crystal; 
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this technique involving making bp = —4/3^ (parameter 
bp defined in Eq.[7|) was discussed in Ref. pL^. In practice, 
however, the effect of the dispersive phase associated with 
the crystal, on the degree of entanglement tends to be 
small. For example, for the KDP source discussed above, 
the value of the Schmidt number, drops from K = 1.065 
to = 1.061, under the effects of chirp compensation. 

In addition to the chirp compensation described above, 
pump chirp can be used as an effective tool to control 
the degree of entanglement in PDC two-photon sources. 
The effect of pump chirp on the JTI may be appreciated 
in Figs. |5|c) and (d). As the pump chirp parameter is 
increased, the correlations apparent in the JTI become 
more pronounced. If the two-photon state is factorable 
in the absence of dispersion, then by adding a controlled 
level of chirp to the pump pulses (prior to the crys- 
tal), it becomes possible to generate two-photon states 
with an arbitrary degree of entanglement, determined by 
the pump chirp /3p. In order to illustrate this behavior. 
Fig. Isja) shows, for a source similar to that assumed for 
Fig. Qb) the Schmidt number, calculated numerically, 
as a function of the pump chirp parameter /3p (where 
we have assumed no dispersion experienced by the sig- 
nal/idler photons). We can see that the Schmidt num- 
ber exhibits a monotonically increasing dependence on 
the pump chirp parameter. Figs, jsjb) through (d) show 
the JTI calculated numerically for three different values 
of the pump chirp, identified along the K vs f3p plot in 
Fig. [sja). Of course, according to Eqns. 



28 



a greater de- 
gree of entanglement as quantified by translates into 
a lower single heralded photon purity, and a lower vis- 
ibility in a two-source HOMI experiment, as demanded 



by the conclusions of Sec. Ill 



The entanglement contributed to the photon pair by 
pump chirp in general manifests itself in the form of mod- 
ified temporal correlations. In the specific case of Fig. |5) 
the strength of temporal correlations increases monoton- 
ically with the pump chirp parameter /3p. These correla- 
tions may be quantified for example through the so-called 
Fedorov ratio [25 , here defined in the temporal domain 



Atc(to) 



(32) 



given in terms of Ate {to) which represents the signal- 
mode temporal width conditioned on a certain idler de- 
tection time to and the unconditional signal-mode tem- 
poral width Ats- Thus, while Ate {to) is defined as the 
full width at half maximum of the function fi{ti) = 
|/(t^,to)P for fixed to selected to maximize the proba- 
bility of emission, Ats is defined as the full width at half 
maximum of the marginal distribution J dts\f{ti^ts)\^. A 
numerical calculation of the Fedorov ratio for the situa- 
tions depicted in Fig.[5|b)-(d) yields the values 1.06, 1.14 
and 3.98 respectively. 



SPECTRAL (TEMPORAL) PROPERTIES OF 
HERALDED SINGLE PHOTONS 



A key application of spontaneous PDC is the genera- 
tion of heralded single photons. Here we consider how 
the continuous variable entanglement properties of the 
photon pair translate into the properties of the heralded 
single photons, in the presence of dispersive elements 
on the paths of the pump, signal and idler. The spec- 
tral (temporal) properties of signal photons heralded by 
an idler detection event may be conveniently character- 
ized in terms of the heralded single photon chronocyclic 
Wigner function (CWF). As shown in Ref. [26], the CWF 
for PDC light characterized by a JSA f{vi^Vs) is given 

by 



Wa{v,t) = J dvog{vo) J div'f {i'o,i'+ y 



X /* ( J^o, - ^ 1 e- 



(33) 



where g{u) represents the idler detection efficiency. The 
CWF is a real- valued function which represents a quasi- 
probability distribution for the emission of a single pho- 
ton with frequency detuning u and at time t] the marginal 
distributions of the CWF yield the corresponding spec- 
tral and temporal single-photon intensity profiles. Here 
we generalize the result in Ref. [26 to the case where the 
three fields involved experience dispersion. We assume 
that the idler detection efficiency has a Gaussian spec- 
tral shape g{u) = exp[— i/^/a^] (centered at = cc^c and 
with bandwidth ag). Carrying out the integrals in Eq.[33| 
with the JSA given in terms of a linear approximation for 
the phase mismatch (see Eq. [4| and using the Gaussian 
approximation for the phasematching function (see text 
before Eq. 15), we obtain 



/ X V4 - T^At^Acj^ 
2.AtA. 
X exp [Tut] , 



r 1 






Aw2_ 


exp 





(34) 



in terms of a mixed-term coefficient F, the temporal du- 
ration At and the spectral width Auj 



r - M 



(35) 



Alo' 



n 



which are expressed, in turn, in terms of the following 
definitions 
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FIG. 5: (a)This plot shows the behavior of the Schmidt number, which quantifies the degree of entanglement, as a function 
of the pump chirp parameter, for the same source as in Fig. [T|a) and (b); this illustrates that pump chirp may be used as a 
powerful tool to control the degree of entanglement. Panels (bj-(d) show the joint temporal intensity, for three different values 
of pump chirp: Pp = Os^, pp = 2.07 x 10~^^s^ and pp = 1.01 x 10~^^s^, identified along the curve in panel (a). 



f2 



rp2 



rp2 



/n~i2 rji2 rTi2 rTi2 \ 

y-'-ssJ-'-ii^R ~ -^si.R-^siJ) 



{Tl. 



,R-'-ii,R ~ \-'-si,R) )' 



(36) 



Here, we have redefined Ta, now referred to as Ta, to 
include the effect of the idler detection efficiency; it is 
expressed as follows 



rp2 rpZ 



t2 ■ 



(37) 



Interestingly, dispersion experienced by the idler pho- 
ton does not have an effect on the spectral (temporal) 
properties of a heralded single photon in the signal mode; 
indeed Eqns. [36] exhibit no dependence on jSi. Note 
also that in the absence of external dispersion, i.e. if 
= Pi = = 0, Tl reduces to = 2T2^^, and 
TI = 0; on the other hand, T3 remains unaffected by the 
presence or absence of dispersion. This in turn implies 
that, according to Eqns. [36} in the absence of dispersion, 
the mixed term F vanishes, in which case we recover the 
expression for the CWF in Ref. [26 . The effect of this 
mixed term, due to dispersion, is that the CWF becomes 
tilted in {cj, t} space in such a way that the frequency 
marginal distribution remains unaffected and the tempo- 
ral marginal distribution is broadened. In order to see 
this explicitly, let us obtain expressions for the marginal 
distributions. The spectral intensity profile of the her- 
alded single photon is given by 



exp 



47rAw2 

in terms of the single-photon spectral width Awm 



1 



(38) 



(39) 



The temporal intensity profile of the heralded single 
photon is given by 



t{t) = j dvWsM 



4 - r^At^Aw^ 



1/2 



exp 



47rAt2 

in terms of the single-photon temporal width A^m 



(40) 



At 



M 



j,2 



(41) 



As expected, while the single-photon spectral width 
Acjm exhibits no dependence on the external disper- 
sion, in general the single-photon temporal width A^m is 
broadened due to dispersion. In order to further under- 
stand this broadening, it is convenient to express the sin- 
gle photon temporal duration as At|^ = Atg + A where 
Atg = 2Tg^g ^ represents the single photon duration for 
no dispersion. We can express A as 



A 

2" 



rp2 (rp2 \2 I rp2 ( rp2 \2 

-'-ii,R\-'-ss,l) -'- ss,R\-'- si,l) 



^-^si,R-^ss 



rp2 rp2 

^ I-'-siJ 



7^2 7^2 
-^ss^R-^iiR 



(42) 
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Re-writing the above expression in terms of the cor- 
relation coefficient S, and using: i)the fact that S is 
constrained by — 1 < S < 1 and ii) the inequality 

-h - 2xy^ > 0, valid for ah x,y with —1<£,<1, 
we can easily show that A > 0, so that Ato represents 
the shortest possible single-photon temporal duration. In 
other words, as expected, the shortest single-photon tem- 
poral duration occurs for no dispersion in any of the three 
optical fields. 

In order to illustrate the characterization of spec- 
tral (temporal) heralded single photon properties. Fig. [g] 
shows a plot of the resulting chronocyclic Wigner func- 
tion for the source parameters assumed in Fig.[l](a), and 
assuming no dispersion experienced by the pump, signal 
and idler modes, no PDC filtering {ctf oo) and ideal 
triggering {dg oo). While the main plot in panel (a) 
shows the numerically-obtained CWF (through numeri- 
cal integration of Eq. 33 ) , the inset shows the CWF plot- 
ted from our analytic expression (Eq. 34). Panels (b) 
and (c) show the CWF marginal distributions, equiva- 
lent to the spectral and temporal single-photon intensity 
profiles. The effects of dispersion on the heralded single- 
photon properties will be discussed in the following sec- 
tion. 
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expressed as f{vi^Vs) = 5{us)I{j^i), the two-photon state 
remains factorable despite the presence of dispersion. In- 
deed, in this case the JSA can then be expressed as: 

f{vi,Vs) = (^(i^s)/(i^i)exp[i(/3p + /3s)i^^] 

X exp[z(/3p + A)^^^]exp[^2/3p^^,^^^]. (43) 

The above expression includes a phase term which is 
not factorable into signal and idler contributions. How- 
ever, the effect of this phase term which depends on the 
product h's'i^ii and which is controlled by the pump chirp 
parameter is suppressed in the limiting case where 
one of the two photons produced is monochromatic; dis- 
persive effects fundamentally require a bandwidth. Note, 
however, that while this implies that an ideal factorable 
two-photon state remains factorable despite an arbitrary 
level of pump chirp, the two photons do individually ex- 
perience temporal broadening due to the phase terms 
proportional to and vf. We can also see from Eq. 43 
that if pump chirp is compensated by an equal magnitude 
of signal and idler chirp, but with the opposite sign, then 
the two photon state becomes unaffected by the presence 
of external dispersion, since in this case phase terms are 
either not present or yield no effect on the two-photon 
state. 

In order to study these effects for a realistic state, i.e. 
one where the idler photon has a small, but non-zero 
bandwidth, let us use our expressions for the two-photon 
state based on the linear dispersion, and Gaussian ap- 
proximations. We will begin this analysis by writing the 
joint spectral amplitude in terms of the adimensional de- 



tunings defined by 



(With/i = 5,z) 



fn{ni,ns) = A^a;exp 



-n? 



2i- 



- n: 



(rp2 rp2 

y ss,R ii,R 



UsUi]. 



)2 



(44) 



FIG. 6: (a) Chronocyclic Wigner function plotted for a source 
identical to that assumed for Fig.[T]^a) and (b), computed nu- 
merically from E g. |33| Inset: plot resulting from the analytic 
expression in Eq. |34| (b) Spectral intensity of heralded single 
photon, in the signal mode, computed as one of the marginal 
distributions of the CWF. (c) Temporal intensity of heralded 
single photon, in the signal mode, computed as one of the 
marginal distributions of the CWF. 



VI. DISPERSION SUPPRESSION EFFECTS 

The factorable, spectrally elongated state leads to 
some interesting dispersion suppression properties. As 
was already discussed in Ref. [l] , in the limit of an "ideal" 
factorable, spectrally elongated state where the signal 
photon is monochromatic, i.e. where the JSA can be 



For the two-photon state to be factorable, the two 
mixed terms (one real and one imaginary) proportional 
to risTii must be vanishingly small. In the absence of ex- 
ternal dispersion, a factorable, asymmetric two-photon 
state can be obtained if: i) an asymmetric group veloc- 
ity matching condition is fulfilled (leading to Ts = 0), 
and ii) for a relatively long crystal, coupled with a rel- 
atively large pump bandwidth, such that \(JTi\ ^ 1. 
Under these conditions, it is straightforward to show 
that the real parts of the T coefficients which define the 
may be expressed as Tf^ ^ 



JSA (see Eq. 



-'- ss,R 



16) 



and 



rp2 



Thus, the correlation 
coefficient becomes S = 2/{^(7Ti) which vanishes for 
the conditions which define this state. The coefficient 
which defines the imaginary mixed term may be written 
as (2(j'^ Pp) / {^(jTi). Thus, a state fulfilling asymmet- 
ric group velocity matching becomes factorable in the 
presence of pump chirp if in addition to \(JTi\ ^ 1, the 
condition \(7Ti\ ^ \cr'^Pp\ is also fulfilled. This last condi- 
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tion tells us that if the adimensional group velocity mis- 
match aTi greatly exceeds the adimensional pump chirp 
the two photon state remains factorable despite 
the presence of pump chirp. Clearly, the larger the co- 
efficient Ti (linearly proportional to the crystal length), 
the more pump chirp can be present while retaining two- 
photon factor ability. Also, note that this "immunity" to 
pump chirp is weakened for larger values of the pump 
bandwidth a. This effect makes factorable and spec- 
trally elongated states, for which the above dispersion 
suppression effect occurs, attractive for practical imple- 
mentations of quantum information processing protocols, 
since in most realistic experimental situations there is 
some chirp present in the pump beam which under typi- 
cal conditions would suppress factorability. 

Recalling that T^^ j = —{/3p-\-/3s) (and similarly for the 



corresponding idler quantity), we can see from Eq. 43 
that if —/3p = /3s = Pi, in addition to the fulfilment 
of the conditions from the previous paragraph, then the 
two-photon state becomes completely immune to exter- 
nal dispersion. In particular, the temporal duration 
of the signal photon in the presence of chirp becomes 
identical to the corresponding duration without exter- 
nal dispersion. In order to see this, it is clear from 



Eq. 

A/A 
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that if Pp = —Ps leading to T, 



2 

ss,I 



- 0, then 
1) (see Eq. l42|. Thus, if 



\crri\ ^ 1 and \cr'^Pp\ <C |crr^| (the conditions discussed in 
the previous paragraph), A ^ and therefore the signal 
temporal duration becomes equal to that in the absence 
of external dispersion. 

In Fig. [7l we illustrate the interplay of pump and sig- 
nal/idler chirp in determining the resulting photon-pair 
temporal properties, for the particular case of a source 
based on a KDP crystal and which fulfils asymmetric 
group velocity matching, similar to that assumed for 
Fig. TJa) and (b). We present plots of the joint tem- 
poral intensity |/(ti,ts)p (panels (a)-(c)), and of the 
chronocyclic Wigner function for the heralded single pho- 
ton (panels (d)-(f)), for a number of different dispersion 
regimes. These plots have been computed numerically 
from Eq. |2] and Eq. [33j without resorting to approxima- 
tions. The insets show the joint temporal intensity cal- 
culated from our analytic expressions (from Eqns. 12 and 
13, where the plots have been centered at ts = ti = 0). 
Note the excellent agreement between the main plots and 
those in the insets; this tells us that the approximation 
used in deriving our analytic expression for the joint tem- 
poral amplitude, i.e. Eq. [T2j involving expressing the 
phasemismatch up to ffist order in frequency detunings, 
is well justified. For these plots we have assumed a fixed 
pump chirp Pp = — 4.77x 10~^^s^ and a varying degree of 
signal and idler chirp (for panels (a) and (d), pg = Pi = 0; 
for panels (b) and (e), ps = Pi = —pp] for panels (c) 
and (f), Ps = Pi = —2Pp). Note that the joint spectral 
intensity |/(i^i,i^s)P exhibits no dependence on external 
dispersion, i.e. it remains unchanged under the effect of 
signal/idler dispersion. 

For this source, the dispersion suppression effect de- 



scribed above may be observed. This effect is clear from 
panels (a)-(c) and from panels (d)-(f). Indeed, the tem- 
poral duration of the signal photon indicated by panel 
(b), represents the minimum which can be obtained (for 
the specific level of pump chirp which has been assumed) , 
and is very close to the temporal duration which would 
be expected in the absence of pump dispersion. This 
can be appreciated from Fig. (sja), where the temporal 
duration AtM of the signal photon is plotted vs the dis- 
persive propagation distance (propagation is assumed to 
occur in fused silica fiber). The dotted blue line indi- 
cates the temporal duration expected without dispersion. 
It is clear that there is a specific value of the propaga- 



tion distance, 



which corresponds to the dispersion 



matching condition ps = Pi = —pp for which the tempo- 
ral duration reaches its minimum value. In Fig. (sjb) we 
show a close-up of the plot shown in panel (a), in the re- 
gion where the minimum occurs. For this source, asym- 
metric group velocity matching leads to = 0, while 
\crTi\ = 1.34 X 10^ so that the condition \crTi\ ^ 1 may be 
considered to be fuffilled. The quantities \cr'^Pp\ and \ari\ 
are of the same order (so that condition \(j'^pp\ <C |crr^| 
is not fulfilled); therefore, for a smaller pump dispersion, 
or a longer crystal, the minimum temporal duration can 
further approach the ideal value obtained without dis- 
persion. 

The effect of chirp on the chronocyclic Wigner func- 
tion may be appreciated from panels (d) and (f) of Fig.[7| 
Specifically, as already discussed in Section [V| chirp in 
general leads to a rotation of this function in chronocyclic 
space in such a way that the spectral width remains con- 
stant, while the temporal width increases. Note that for 
the dispersion matching propagation distance, for which 
Pp = —Ps = —Pi, the structure of the CWF and the JTI 
without dispersion is essentially recovered, despite the 
presence of external dispersion (i.e. compare Fig. [Tj^b) 
with Fig.jljb), and Fig. [T^e) with Fig.jeja)). 



VII. EFFECTS OF DISPERSIVE 
PROPAGATION ON THE TWO-PHOTON STATE 



In this section we explore the effect on the spectral 
characteristics of PDC two-photon states of propagation 
through a dispersive medium. We will consider the spe- 
cific case where the signal and idler photons both prop- 
agate through the same dispersive medium, such as a 
fiber, assumed to be non-birefringent so that the two 
orthogonally-polarized generated photons experience the 
same dispersion (these results could be easily generalized 
to differing dispersion for the signal and idler photons, 
since our treatment in previous sections is valid for both 
balanced and unbalanced propagation). 

The two photon state, after propagation through dis- 
tance z (e.g. in a fiber) in a medium with quadratic 
chirp (i.e. ignoring higher order dispersive terms), can 
be written as 
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FIG. 7: Panels (a)-(c): Joint temporal intensity for three different levels of chirp [/3 = 0,/3 = —(3p and (3 = — 2/3p, corresponding 
to the following lengths of propagation through fused silica fiber: L = Om, L = 2.64m and L = 5.28m] experienced by the signal 
and idler photons, and a pump chirp level of f3p = —4.77 x 10~^^s^. While these plots were computed numerically through 
Eq. [2] the insets show the corresponding plots computed from our analytical expression, see Eq. |12| (d )-(f ): Chronocyclic 
Wigner function corresponding to panels (a)-(c). While these plots were computed numerically through Eq.[33] the insets show 
the corresponding plots computed from our analytical expression, see Eq. |34[ 



l^^iz)) = |vac)+7? J dvs j dvif{vi,Vs]z)\vs)\vi) (45) 
in terms of the corresponding JSA 



(46) 



Here, /(i^i, I's]^) denotes the JSA for z = i.e. for no 
dispersion experienced by the signal and idler photons 
{B = 0), as given by Eq. 15 B is the GVD parameter 
defined as B = hc" /2 where represents the wavenumber 
which characterizes the dispersive medium. Note that 
/(i^i, Vs] 0) includes any dispersive phases associated with 
pump chirp. 

While the degree of entanglement remains constant 
during propagation of the signal/idler photons through a 
lossless medium, we will show that the degree of modulus- 
only entanglement, i.e. ignoring any phase entanglement, 
when calculated in the time domain can vary drastically 



during propagation. What this means is that entangle- 
ment can be considered to reside to varying degrees in the 
modulus and in the phase of the joint temporal ampli- 
tude depending on the propagation distance z. We refer 
to this phenomenon as spectral (temporal) entanglement 
migration, referring to the observation that entanglement 
can "migrate" between the modulus and the phase of the 
JTA; there is a corresponding effect which occurs in the 
spatial domain for free-space propagation [27] . 

In our analysis we first restrict attention to the mod- 
ulus of the joint amplitude, i.e. we neglect any phase 
entanglement. In cases where there is in fact no phase 
entanglement, then of course the degree of entanglement 
of the reduced state considered here coincides with the 
degree of entanglement of the actual physical state. Also 
note that the joint spectrum obtained through a phase- 
insensitive spectrographic measurement such as those re- 
ported in Refs. [28l |29] corresponds to the reduced state 
considered here. By carrying out a Schmidt decompo- 
sition of the reduced sate and computing the resulting 
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FIG. 8: (Color online) (a) Temporal duration of heralded 
single photon in the signal mode, as a function of the chirp 
parameter which characterizes the propagation of signal and 
idler photon pairs. While the dashed black line was computed 
through a numerical calculation, the re d lin e was computed 
through our analytic expression (see Eq. 41 ). In this case the 
pump is assumed to be chirped with f3p = —4.77 x 10~^^s^. 
Note that the temporal duration as a function of chirp pa- 
rameter exhibits a minimum at the propagation distance cor- 
responding to /3 = —Pp- The dotted blue line represents the 
temporal duration for no external dispersion, (b) Close-up of 
the previous panel in the region of the minimum. 



Schmidt number, for different values of the propagation 
distance it becomes possible to directly monitor the de- 
gree of modulus-only entanglement during propagation of 
the photon pair. We refer to the modulus-only Schmidt 
number calculated in the spectral domain as Km,s and 
to the modulus-only Schmidt number in the temporal 
domain as Km,T- Note that while the Schmidt num- 
ber K calculated in the temporal domain must equal the 
Schmidt number calculated in the spectral domain, this 
is no longer true for the modulus-only Schmidt numbers 
where, in general, Km,T 7^ Km.s- We also note that 
while Km^s is constant with respect to the propagation 
distance z (since the JSI is insensitive to dispersion), in 
general Km^T is a function of z. 

While carrying out an analytic Schmidt decomposi- 
tion on the full two-photon state without resorting to 
approximations is challenging, this becomes possible for 
two-photon states for which the JSA can be written in 
terms of a real-valued Gaussian function of the form 
exp[— (Ax^ + By^ + Cxy)] where x^y are variables and 
A, 5, C parameters. In this case, the Schmidt eigen- 



functions are given by Gauss-Hermite functions (see 
Refs. [22l |30]), while the eigenvalues are given as = 
(1 — /i^)/i^^ characterized by a parameter /i. Calculation 
of the Schmidt decomposition in the temporal domain 
yields the following expression for this parameter 



= ' (47) 

in terms of which the reduced Schmidt number (in the 
temporal domain) is given by 



2 • 



(48) 



In order to illustrate this entanglement migration ef- 
fect, we will consider as specific example the factorable 
(in the absence of external dispersion) states obtained 
through asymmetric GVM. 
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FIG. 9: (Color online) This plot shows, as a function of 
the propagation distance the Schmidt number K and the 
reduced temporal Schmidt number Km,T] we present both 
numerical (dashed black line) and analytic (continuous red 
line) curves for these functions. The minimum allowed value 
K — 1 IS indicated with a dotted blue line. While K remains 
constant during propagation, Km,T exhibits a rich structure. 
From these plots we may infer that entanglement migrates 
from the modulus of the joint temporal amplitude (at z = zi), 
to its phase (at z = Z2), back to the modulus (at 2; = zs), and 
slowly to the phase once more (for large z). 

Fig. [9] shows, for a factorable, spectrally elongated 
source (with the same parameters as in Fig. [ija)), the 
dependence of the reduced Schmidt number Kjn,T^ cal- 
culated through Eq. |48) on the propagation distance z 
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(through a fused sihca fiber). Of course, the dependence 
of the "total" Schmidt number K on propagation dis- 
tance z is trivial: it remains constant. Here, we have 
assumed a value for the quadratic pump dispersion chirp 
of Pp = —4.77 X 10~^^5^ (note that for vanishing pump 
chirp, there is no entanglement, and therefore Km,T re- 
mains constant at unity). We observe that there is a spe- 
cific propagation distance, Zmin^ for which the reduced 
Schmidt number Km,T essentially reaches the minimum 
value possible, i.e. unity, which corresponds to the dis- 
persion matching length for which /3s = /3i = —/3p iden- 
tified in the previous section. In the figure we have pre- 
sented both: i) a plot of Km,T derived from our ana- 



lytic expression (continuous red line; see Eq. 48), and 
ii) a numerical version where we have not resorted to 
approximations (dashed black line). For the analytic cal- 
culation, we obtain the "total" K through Eq. |30) with 
K = 1/V. The plots obtained analytically and numer- 
ically exhibit the same general features for Km,T{z): a 
fast drop from z = to z = Zmin (corresponding to 
P = — /^p), reaching a value approaching unity, a fast rise 
to a level similar to that at z = 0, and a gradual drop 
for larger propagation distances. Thus, in the temporal 
domain, as the photon pair propagates through the fiber, 
entanglement migrates from the modulus to the phase, 
back to the modulus, and finally gradually to the phase 
once more. This behavior is, of course, also apparent 
from Fig. [t] panels (a)-(c), where the first and last pan- 
els show strong temporal (modulus) correlations, while 
the dispersion-matched case shows essentially no tempo- 
ral (modulus) correlations. Note that a similar behavior 
to that observed in Fig. |9] for Km,T would be observed in 
terms of the Fedorov ratio (see Eq. 32). 

Let us note that for the dispersion-matching propaga- 
tion distance in the fiber, z = Zmim Km,T ~ 1 implies 
that in the time domain entanglement resides entirely 
in the phase. Furthermore, as discussed above for this 
particular source, when viewed in the spectral domain, 
entanglement resides entirely in the phase for any prop- 
agation distance z. Thus, remarkably, for this particular 
source at z = Zmin entanglement resides entirely in the 
phase in both domains. We speculate that applications of 
this form of phase-only entanglement, in both the spec- 
tral and temporal domains, may result from the fact that 
it is less accessible to measurements than more typical 
modulus entanglement. 



pansion of the phasemismatch, as well as a corresponding 
expression where we also resort to the Gaussian approx- 
imation for the phasematching function. While the first 
one yields two-photon characteristics remarkably simi- 
lar to those obtained through a numerical calculation 
which does not resort to approximations, the second one 
permits analytic expressions for the following: i) Hong- 
Ou-Mandel interferogram, involving two heralded single 
photon sources, ii) Chronocyclic Wigner function for her- 
alded single photons, which fully characterizes the spec- 
tral (temporal) properties of the single photons and iii) 
Schmidt decomposition of the modulus of the joint tem- 
poral amplitude. Exploiting this description, we find that 
pump chirp can represent an extremely useful tool for 
controlling the degree of photon-pair entanglement. In- 
deed, we find that if the two-photon state is factorable 
in the absence of pump chirp, then the resulting Schmidt 
number (which quantifies the entanglement present) can 
be continuously adjusted, in principle, from K = 1 to 
any desired value through the pump chirp parameter. We 
likewise predict a dispersion suppression effect which oc- 
curs for two-photon states designed to be factorable and 
spectrally-elongated: the temporal duration of the her- 
alded single photon in the presence of pump chirp can 
approach the temporal duration which would have been 
observed without pump chirp, if a specific condition is 
fulfilled by the pump and PDC chirp parameters. 

We explicitly show the relationship between: i) the 
photon pair Schmidt number, ii) the Hong-Ou-Mandel 
visibility in a two-crystal arrangement and iii) the her- 
alded single photon purity. This gives us a method to ex- 
perimentally characterize both the degree of photon-pair 
entanglement and the degree of heralded single photon 
purity. We show that the spectral (temporal) entangle- 
ment present in a two photon state can migrate between 
the modulus and the phase of the joint temporal ampli- 
tude. We present details of a specific source for which 
entanglement is expected to reside entirely in the phase 
in both the spectral and temporal domains. We believe 
that the analysis presented in this paper will be useful 
for the design and implementation of specific photon-pair 
sources tailored for quantum information processing ap- 
plications. 
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In this paper we have presented a study of the effects 
of external dispersion, i.e. introduced by optical elements 
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sion. We have derived expressions for the joint amplitude 
of the photon pair state, both in the spectral and tempo- 
ral domains. In particular, we present an expression for 
the joint temporal amplitude resorting to a first-order ex- 
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APPENDIX A 

The analytic expression for the joint amplitude (see 
Eq. 12) may be written in terms of a function of the 
form 
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G{x-\-xo) 

Z{G,xo;x) = ^ J dte-'" 

G{x — xo) 

= \ {erf[G(x + xo)] - erf[G(x - xq)]} (A1) 

If we write parameter G in terms of its modulus and 
phase, i.e. G = |G'|exp(z^), this function may be shown 
to converge for — 7r/4 < < tt/A. 

In the hmit where |Gxo| oo, function Z{G,xo;x) 
behaves as a "top hat" function 

hm |Z(G, xo; x)| = rect(— xo, xo; x) (A2) 

\Gxo\^oo 



where 

rect(x;a,6) = |Q^ otherwise ^^^^ 

In the hmit where \Gxo\ 0, function Z{G^xo;x) 
behaves as a Gaussian function. 



hm Z{GjXo;x) = —^Gxoex.-p{—G'^x^) (A4) 
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